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Abstract. The Calculus of Wrapped Compartments is a framework 
based on stochastic multiset rewriting in a compartmentalised setting 
originally developed for the modelling and analysis of biological interac¬ 
tions. In this paper, we propose to use this calculus for the description 
of ecological systems and we provide the modelling guidelines to encode 
within the calculus some of the main interactions leading ecosystems 
evolution. As a case study, we model the distribution of height of Croton 
wagneri , a shrub constituting the endemic predominant species of the 
dry ecosystem in southern Ecuador. In particular, we consider the plant 
at different altitude gradients (i.e. at different temperature conditions), 
to study how it adapts under the effects of global climate change. 
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1 Introduction 

Answers to ecological questions could rarely be formulated as general laws: ecol¬ 
ogists deal with in situ methods and experiments which cannot be controlled in 
a precise way since the phenomena observed operate on much larger scales (in 
time and space) than man can effectively study. Actually, to carry on ecological 
analyses, there is the need of a “macroscope”! 

Theoretical and Computational Ecology, the scientific disciplines devoted to 
the study of ecological systems using theoretical methodologies together with 
empirical data, could be considered as a fundamental component of such a 
macroscope. Within these disciplines, quantitative analysis, conceptual descrip¬ 
tion techniques, mathematical models, and computational simulations are used 
to understand the fundamental biological conditions and processes that affect 
populations dynamics (given the underlying assumption that phenomena ob¬ 
servable across species and ecological environments are generated by common, 
mechanistic processes) [57i . 

Ecological models can be deterministic or stochastic [18]. Given an initial 
system, deterministic simulations always evolve in the same way, producing a 

* This work has been partially sponsored by the BioBITs Project (Converging Technologies 2007 , 
area: Biotechnology-ICT), Regione Piemonte. Angelo Troina has visited the Instituto de Ecologfa 
at UTPL (Loja, Ecuador) within the Prometeo program founded by SENESCYT. A preliminary 
version of this paper has been presented in 1601 . 
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unique output [65] • Deterministic methods give a picture of the average, expected 
behaviour of a system, but do not incorporate random fluctuations. On the other 
hand, stochastic models allow to describe the random perturbations that may 
affect natural living systems, in particular when considering small populations 
evolving at slow interactions. Actually, while deterministic models are approxi¬ 
mations of the real systems they describe, stochastic models, at the price of an 
higher computational cost, can describe exact scenarios. 

A model in the Calculus of Wrapped Compartments (CWC for short) con¬ 
sists of a term, representing a (biological or ecological) system and a set of 
rewrite rules which model the transformations determining the system’s evolu¬ 
tion 1271231 - Terms are defined from a set of atomic elements via an operator of 
compartment construction. Each compartment is labelled with a nominal type 
which identifies the set of rewrite rules that may be applied into it. The CWC 
framework is based on a stochastic semantics and models an exact scenario able 
to capture the stochastic fluctuations that can arise in the system. 

The calculus has been extensively used to model real biological scenarios, in 
particular related to the AM-symbiosis pugg An hybrid semantics for CWC, 
combining stochastic transitions with deterministic steps, modelled by Ordinary 
Differential Equations, has been proposed in [25126] . 

While the calculus has been originally developed to deal with biomolecular in¬ 
teractions and cellular communications, it appears to be particularly well suited 
also to model and analyse interactions in ecology. In particular, we present in 
this paper some modelling guidelines to describe, within CWC, some of the main 
common features and models used to represent ecological interactions and popu¬ 
lation dynamics. A few generalising examples illustrate the abstract effectiveness 
of the application of CWC to ecological modelling. 

As a real case study, we model the distribution of height of Croton wagneri , 
a shrub in the dry ecosystem of southern Ecuador, and investigate how it could 
adapt to global climate change. 


1.1 Motivation and Methodology 

At the beginning of the second half of the twentieth century, when ecology was 
still a young science and mathematical models for ecological systems were in their 
infancy, Elton i.35], acknowledging the influence of Lotka [44] and Volterra in], 
wrote: “Being mathematicians, they did not attempt to contemplate a whole 
food-chain with all the complications of five stages. They took two: a predator 
and its prey”. Nowadays, in the era of computational ecological modelling, de¬ 
terministic systems based on ordinary differential equations for two variables, 
or even a whole food chain, appear like simple idealisations quite distant from 
the real complexity of nature. Predator-prey interactions are now considered as 

3 Arbuscular Mycorrhiza (AM) is a class of fungi constituting a vital mutualistic 
interaction for terrestrial ecosystems. More than 48% of land plants actually rely on 
mycorrhizal relationships to get inorganic compounds, trace elements, and resistance 
to several kinds of pathogens. 










Ecological Modelling with the Calculus of Wrapped Compartments 


3 


“consumer resource” interactions embedded within the large ecological networks 
that underlie biodiversity. Lotka-Volterra equations and their many descendants 
assume that individuals within a system are well mixed and interact at mean 
population abundances. They are mean-field equations that use the mass-action 
law to describe the dynamics of interacting populations, and ignore both the 
scale of individual interactions and their spatial distribution. However, because 
ecological systems are typically nonlinear, they often cannot be solved analyti¬ 
cally, and, in order to obtain sensible results, nonlinear, stochastic computational 
techniques must be used. 

The formal framework to be used as the modelling core of this project should 
thus be able to manage several features which are typical of ecological sys¬ 
tems. Namely, complex ecological systems are multilevel, they follow non linear, 
stochastic dynamics and involve a distributed spatial organisation. 

Multilevel Modelling. The role of the computational methodology used to 
model and simulate ecological systems is to address questions on the relation¬ 
ship between systems dynamics at different temporal, spatial, and organisational 
(or structural) scales. In particular, it is important to address the variability at 
small, local scales and its effects on the dynamics of the aggregated quantities 
measured at large, global scales [54] . 

Stochastic Modelling. Ecological models can be deterministic or stochas¬ 
tic |18j . Given an initial system, deterministic simulations always evolve in the 
same way, producing a unique output [55] . Deterministic methods give a pic¬ 
ture of the average, expected behaviour of a system, but do not incorporate 
random fluctuations. On the other hand, stochastic models allow to describe 
the random perturbations that may affect natural living systems, in particu¬ 
lar when considering small populations evolving at slow interactions. Actually, 
while deterministic models are approximations of the real systems they describe, 
stochastic models, at the price of an higher computational cost, can describe ex¬ 
act scenarios. Stochastic models, such as interacting particle systems, can also 
help us examine new approaches for scaling up individual-based dynamics]^] 

Spatial Modelling. The impact of space tends to make the population dynam¬ 
ics significantly more complicated compared with its non-spatial counterpart 
and to bring new and bigger challenges to simulations. Formal models dealing 
explicitly with spatial coordinates are able to depict more precise localities and 
ecological niches , describing, for example, how organisms or populations respond 
to the distribution of resources and competitors H51 . 

4 Note that the impact of stochastic factors and the corresponding level of a system’s 
uncertainty are much higher in Ecology than in other natural sciences. Common 
sources of uncertainty are, e.g., the poor accuracy of ecological data and their tran¬ 
sient nature. Noise that is inevitably present in ecosystems can significantly change 
the properties of an ecological model and this fundamental uncertainty affects the 
accuracy of ecological data [55l . 
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The Calculus of Wrapped Compartments. While the Calculus of Wrapped 
Compartments has been originally developed to deal with biomolecular inter¬ 
actions and cellular communications, it appears to be particularly well suited 
also to model and analyse interactions in ecology. The Calculus of Wrapped 
Compartments satisfies the main requirements addressed in the previous sec¬ 
tion. Namely, CWC is able to model and simulate: (i) multilevel systems, (ii) 
stochastic dynamics, (iii) explicit spatial systems. 

The compartment operator of the calculus can be used to describe the topo¬ 
logical organisation of a systems. It also allows to deal with multilevel systems 
by defining different set of rules for different compartments, reflecting the inter¬ 
actions taking place at the different levels of the system. 

The evolution of a system described in CWC follows a stochastic simulation 
model defined by incorporating a collision-based framework along the lines of 
the one presented by Gillespie in [38], which is, de facto, the standard way to 
model quantitative aspects of biological systems. The basic idea of Gillespie’s 
algorithm is that a rate is associated with each considered reaction. This rate 
is used as the parameter of an exponential probability distribution modelling 
the time needed for the reaction to take place. In the standard approach, the 
reaction propensity is obtained by multiplying the rate of the reaction by the 
number of possible combinations of reactants in the compartment in which the 
reaction takes place, modelling the law of mass action. 

A spatial extension of CWC has been proposed in mi. incorporating a two- 
dimensional spatial description of the elements in the system through axial coor¬ 
dinates and special rules for the movement of system components in space. The 
spatial extension of the calculus can be generalised to deal with spaces defined 
in more than two dimensions. 


Summary In Section [2] we recall the syntax and semantics of CWC. In Section [3] 
we present some of the characteristic processes leading ecosystems evolution and 
show how to encode them in CWC. In Section [4] we model the distribution of 
Croton wagneri, a shrub in the dry ecosystem of southern Ecuador, and investi¬ 
gate how it could adapt to global climate change. Finally, in Section [5] we draw 
our conclusions and survey some related work. 


2 The Calculus of Wrapped Compartments 

The Calculus of Wrapped Compartments (CWC) (see [271251201 1 is based on a 
nested structure of compartments delimited by wraps with specific proprieties. 


2.1 Term Syntax 

Let A be a set of atomic elements (atoms for short), ranged over by a, b, ..., and 
£ a set of compartment types represented as labels ranged over by i, (!, l\, ... 
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Definition 1 (CWC terms). A CWC term is a multiset t of simple terms t 
defined by the following grammar: 

t ::= a | (aj t'Y 

A simple term is either an atom or a compartment consisting of a wrap (rep¬ 
resented by the multiset of atoms a), a content (represented by the term t') and 
a type (represented by the label £). Multisets are identified modulo permutations 
of their elements. The notation n * t denotes n occurrences of the simple term t. 
We denote an empty term with •. 

In applications to ecology, atoms can be used to describe the individuals of 
different species and compartments can be used to distinguish different ecosys¬ 
tems, habitats or ecological niches. Compartment wraps can be used to model 
geographical boundaries or abiotic components (like radiations, climate, atmo¬ 
spheric or soil conditions, etc.). In evolutionary ecology, individuals can also be 
described as compartments, showing characteristic features of their phenotype 
in the wrap and keeping their genotype (or particular alleles of interest) in the 
compartment content. 

An example of CWC term is 20*a 12 * b (c d J 6 * e 4 * f) e representing a 
multiset (denoted by listing its elements separated by a space) consisting of 20 
occurrences of a, 12 occurrence of b (e.g. 32 individuals of two different species) 
and an Ctype compartment (c d J 6* e 4* f) e which, in turn, consists of a wrap 
(a boundary) with two atoms c and d (e.g. two abiotic factors) on its surface, 
and containing 6 occurrences of the atom e and 4 occurrences of the atom / 
(e.g. 10 individuals of two other species). Compartments can be nested as in the 
term (a b cj (d ej f) e g h) e . 


2.2 Rewrite Rules 


System transformations are defined by rewrite rules, defined by resorting to 
CWC terms that may contain variables. 


Definition 2 (Patterns and Open terms). Simple patterns P and simple 
open terms O are given by the following grammar: 


P 

O 

q 


a {a x\P Xy 
a 0 q\Oy | X 

a x 


where a is a multiset of atoms, P is a pattern (i.e., a, possibly empty, multiset 
of simple patterns), x is a wrap variable (can be instantiated by a multiset of 
atoms), X is a content variable (can be instantiated by a CWC term), q is a 
multiset of atoms and wrap variables and O is an open term (i.e., a, possibly 
empty, multiset of simple open terms). 


We will use patterns as the l.h.s. components of a rewrite rule and open terms 
as the r.h.s. components of a rewrite rule. Patterns are intended to match, via 
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substitution of variables, with ground terms (containing no variables). Note that 
we force exactly one variable to occur in each compartment content and wrap of 
our patterns. This prevents ambiguities in the instantiations needed to match a 
given compartment]^] 

Definition 3 (Rewrite rules). A rewrite rule is a triple (£, P, O), denoted by 
t : P i —> O, where the pattern P and the open term O are such that the variables 
occurring in O are a subset of the variables occurring in P. 

The rewrite rule £ : P e->• O can be applied to any compartment of type £ 
with P in its content (that will be rewritten with O). Namely, the application 
of t : P K > O to term t is performed in the following way: 

1. find in t (if it exists) a compartment of type £ with content t' and a substi¬ 
tution a of variables by ground terms such that t' = a(P X)0 

2. replace in t the subterm t ' with <j{0 X). 

For instance, the rewrite rule £ : a b i—>■ c means that in compartments of 
type £ an occurrence of a b can be replaced by c. We write t H > t' to denote a 
reduction obtained by applying a rewrite rule to t resulting to t'. 

While a rewrite rule does not change the label £ of the compartment where 
it is applied, it may change the labels of the compartments occurring in its 
content. For instance, the rewrite rule £ : (a £ J X)* 1 ► (a x\ XY 2 means that, 

if contained in a compartment of type £, a compartment of type l\ containing 
an a on its wrap can be changed to type £n- 

CWC Models. For uniformity reasons we assume that the whole system is al¬ 
ways represented by a term consisting of a single (top level) compartment with 
distinguished label T and empty wrap, i.e., any system is represented by a term 
of the shape (*J f) T , which, for simplicity, will be written as t. Note that while 
an infinite set of terms and rewrite rules can be defined from the syntactic defi¬ 
nitions in this section, a CWC model consists of an initial system (*J t) T and a 
finite set of rewrite rules 1Z. 

2.3 Stochastic Simulation 

A stochastic simulation model for ecological systems can be defined by incor¬ 
porating a collision-based framework along the lines of the one presented by 
Gillespie in [3S], which is, de facto , the standard way to model quantitative as¬ 
pects of biological systems. The basic idea of Gillespie’s algorithm is that a rate 
is associated with each considered reaction which is used as the parameter of an 
exponential probability distribution modelling the time needed for the reaction 

5 The linearity condition, in biological terms, corresponds to excluding that a trans¬ 
formation can depend on the presence of two (or more) identical (and generic) com¬ 
ponents in different compartments (see also |52j). 

6 The implicit (distinguished) variable X matches with all the remaining part of the 
compartment content. 
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to take place. In the standard approach the reaction propensity is obtained by 
multiplying the rate of the reaction by the number of possible combinations of 
reactants in the compartment in which the reaction takes place, modelling the 
law of mass action. 

Stochastic rewrite rules are thus enriched with a rate k (notation £ : P i —> 


O). Evaluating the propensity of the stochastic rewrite rule R = l : a b i—>• c 
within the term t = a a a b b, contained in the compartment u = (*J tY, we 
must consider the number of the possible combinations of reactants of the form 
a b in t. Since each occurrence of a can react with each occurrence of &, this 
number is 3 • 2, and the propensity of R within u is k ■ 6. A detailed method to 
compute the number of combinations of reactants can be found in m- 

The stochastic simulation algorithm produces essentially a Continuous Time 
Markov Chain (CTMC). Given a term t, a set 1Z of rewrite rules, a global time 
6 and all the reductions e\,... , eu applicable in all the different compartments 
of t. with propensities n,. .. ,rjn, Gillespie’s “direct method” determines: 


— The exponential probability distribution (with parameter r = Y^iLi r i ) of 
the time r after which the next reduction will occur; 

— The probability ?y/r that the reduction occurring at time S + r will be 


The CWC simulator [2j is a tool under development at the Computer Sci¬ 
ence Department of the Turin University, based on Gillespie’s direct method 
algorithm |3SJ. It treats CWC models with different rating semantics (law of 
mass action, Michaelis-Menten kinetics, Hill equation) and it can run indepen¬ 
dent stochastic simulations over CWC models, featuring deep parallel optimiza¬ 
tions for multi-core platforms on the top of FastFlow [5]. It also performs online 
analysis by a modular statistical framework m- 


3 Modelling Ecological Systems in CWC 

Computational Ecology is a field devoted to the quantitative description and 
analysis of ecological systems using empirical data, mathematical models (in¬ 
cluding statistical models), and computational technology. While the different 
components of this interdisciplinary field of research are not new, there is a 
new emphasis on the integrated treatment of the area. This emphasis is am¬ 
plified by the expansion of our local, national, and international computational 
infrastructure, coupled with the heightened social awareness of ecological and 
environmental issues and its effects on research funding. 

We advocate a convergence between computer and life sciences. This emerg¬ 
ing paradigm moves to a system level understanding of life, where unpredictable, 
complex behaviour show up. We claim that computer science will greatly con¬ 
tribute to a better understanding of the behaviour of ecological systems. We 
plan to develop models, languages and tools for describing, analysing and imple¬ 
menting in silico ecological systems, as an additional contribution of Information 

' Reductions are applied in a sequential way, one at each step. 
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Technology to those typical research areas in current Computational Ecology, 
such as (i) storing, organising and retrieving large amounts of ecological data or 
(ii) visual modelling techniques for scientific visualisation of multi-dimensional, 
computer-generated scenes that can be used to express empirical data. 

More in detail, we use our formal framework for modelling and studying the 
behaviour of living systems. Our starting point is that ecological systems are 
conveniently described as entities that change their state because of the occur¬ 
rence of biotic and abiotic interactions, giving rise to some observable behaviour. 
We thus adhere to the view of living systems as biological computing units. 

In this section we present some of the characteristic features leading the 
evolution of ecological systems, and we show how to encode them within CWC. 


3.1 Population Dynamics 

Models of population dynamics describe the changes in the size and composition 
of populations. 

The exponential growth model is a common mathematical model for pop¬ 
ulation dynamics, where, using r to represent the pro-capita growth rate of a 
population of size N, the change of the population is proportional to the size of 
the already existing population: 



CWC Modelling 1 (Exponential Growth Model) We can encode within 
CWC the exponential growth model with rate r using a stochastic rewrite rule 
describing a reproduction event for a single individual at the given rate. Namely, 
given a population of species a living in an environment modelled by a com¬ 
partment with label t, the following CWC rule encodes the exponential growth 
model: 

i : a i——f a a 

Counting the number of possible reactants, the growth rate of the overall popu¬ 
lation is automatically obtained by the stochastic semantics underlying CWC. 

A metapopulatior^is a group of populations of the same species distributed in 
different patched and interacting at some level. Thus, a metapopulation consists 
of several distinct populations and areas of suitable habitat. 

Individual populations may tend to reach extinction as a consequence of 
demographic stochasticity (fluctuations in population size due to random demo¬ 
graphic events); the smaller the population, the more prone it is to extinction. A 
metapopulation, as a whole, is often more stable: immigrants from one popula¬ 
tion (experiencing, e.g., a population boom) are likely to re-colonize the patches 


The term metapopulation was coined by Richard Levins in 1970. In Levins’ own 
words, it consists of “a population of populations” [42] , 

9 A patch is a relatively homogeneous area differing from its surroundings. 
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left open by the extinction of other populations. Also, by the rescue effect , in¬ 
dividuals of more dense populations may emigrate towards small populations, 
rescuing them from extinction. 

Populations are affected by births and deaths, by immigrations and emi¬ 
grations (BIDE model l 23]). The number of individuals at time t + 1 is given 
by: 

iVt+i = Nt + B + I — D — E 

where N t is the number of individuals at time t and, between time t and t + 1, 
B is the number of births, / is the number of immigrations, D is the number 
of deaths and E is the number of emigrations. Conditions triggering migration 
could be: climate, food availability or mating [33| . 


CWC Modelling 2 (BIDE model) We can encode within CWC the BIDE 
model for a compartment of type I using stochastic rewrite rules describing the 
given events with their respective rates r, i, d, e: 


T : a (x\XY^(x\a X) e 
i : a • 

T : (®J a XY h-4 a (a:J X) 1 


(birth) 

(immigration) 

(death) 

(emigration) 


Starting from a population of N t individuals at time t, the number N t+ \ of indi¬ 
viduals at time t +1 is computed by successive simulation steps of the stochastic 
algorithm. The race conditions computed according to the propensities of the 
given rules assure that all of the BIDE events are correctly taken into account. 


Example 1. Immigration and extinction are key components of island biogeogra¬ 
phy. We model a metapopulation of species a in a context of 5 different patches: 
4 of which are relatively close, e.g. different ecological regions within a small 
continent, the last one is far away and difficult to reach, e.g. an island. The 
continental patches are modelled as CWC compartments of type I c , the island 
is modelled as a compartment of type £,. Births, deaths and migrations in the 
continental patches are modelled by the following CWC rules: 

n 0.005 n 0.005 

i c : a i —> a a i c : a \—> • 

T : (ccj a X) l ° ^ a (sj X) 1 * T : a (a:J X)^ (sj a X )*> 

These rates are drawn considering days as time unites and an average of life 
expectancy and reproduction time for the individuals of the species a of 200 days 
(oho!;) - For t ^ ie m °delling of real case studies, these rates could be estimated from 
data collected in situ by tagging individuals (^] In this model, when an individual 
emigrates from its previous patch it moves to the top-level compartment from 
where it may reach one of the close continental patches (might also be the old 


10 


In the remaining examples we will omit a detailed time description. 
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one) or start a journey through the sea (modelled as a rewrite rule putting the 
individual on the wrapping of the island compartment): 


T : a (xj X) 


L °- 2 , 


(x aj X) ei 


Crossing the ocean is a long and difficult task and individuals trying it will 
probably die during the cruise; the luckiest ones, however, might actually reach 
the island, where they could eventually benefit of a better life expectancy for 
them and their descendants: 


T: (x a\Xy< °h^ 3 {x]X) t * T: {xa\Xf 

„ 0.007 „ 0.003 

Ki : a i —> a a : a i—► • 


0.0005 


(x J a xy 


Considering the initial system modelled by the CWC term: 

t = (,j30*a) 4 (*j30*a) 4 (*j30*a) 4 (*j30>i=a) 4 (•]•)* 

we can simulate the possible evolutions of the overall diffusion of individuals 
of species a in the different patches. Notice that, on average, one over °j^ 0 3 0 3 5 
individuals that try the ocean journey, actually reach the island. In Figure [l] we 
show the result of a simulation plotting the number of individuals in the different 
patches in a time range of approximatively 10 years. Note how, in the final part 
of the simulation, empty patches get recolonised. In this particular simulation, 
also, an exponential growth begins after the colonisation of the island. The full 
CWC model describing this example can be found at: http://www.di.unito. 
it/~troina/cmc!3/metapopulation.cwc 



Fig. 1. Metapopulation dynamics. 
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In ecology, using r to represent the pro-capita growth rate of a population and 
K the carrying capacity of the hosting environment j^Jr/AT selection theory [56] 
describes a selective pressure driving populations evolution through the logistic 
model [70] : 


dN 

dt 


= r ■ N ■ 



N 

I< 


where N represents the number of individuals in the population. 


CWC Modelling 3 (Logistic Model) The logistic model with growth rate r 
and carrying capacity K, for an environment modelled by a compaHment with 
label £, can be encoded within CWC using two stochastic rewrite rules describ¬ 
ing (i) a reproduction event for a single individual at the given rate and (ii) a 
death event modelled by a fight between two individuals at a rate that is inversely 
proportional to the carrying capacity: 


£ 

£ 


r 



a a 


: a a 


a 


If N is the number of individuals of species a, the number of possible reactants 
for the first rule is N and the number of possible reactants for the second rule 
is, in the exact stochastic model, (^) = n '(n-i) , j e num }) er 0 f distinct 
pairs of individuals of species a. Multiplying this values by the respective rates 
we get the propensities of the two rules and can compute the value of N when 
the equilibrium is reached (i.e., when the propensities of the two rules are equal): 
r • N = ^ w h en jV = 0 or N = K. 

For a given species, this model allows to describe different growth rates and 
carrying capacities in different ecological regions. Identifying a CWC compart¬ 
ment type (through its label) with an ecological region, we can define rules 
describing the growth rate and carrying capacity for each region of interest. 

Species showing a high growth rate are selected by the r factor, they usually 
exploit low-crowded environments and produce many offspring, each of which has 
a relatively low probability of surviving to adulthood. By contrast, A"-selected 
species adapt to densities close to the carrying capacity, tend to strongly compete 
in high-crowded environments and produce fewer offspring, each of which has a 
relatively high probability of surviving to adulthood. 

Example 2. There is little, or no advantage at all, in evolving traits that permit 
successful competition with other organisms in an environment that is very likely 
to change rapidly, often in disruptive ways. Unstable environments thus favour 
species that reproduce quickly (r-selected species). Characteristic traits of r- 
selected species include: high fecundity, small body, early reproduction and short 
generation time. Stable environments, by contrast, favour the ability to compete 


ii 


I.e., the population size at equilibrium. 
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successfully for limited resources (A'-selected species). Characteristic traits of K- 
selected species include: large body size, long life expectancy, production of fewer 
offspring (usually requiring extensive parental care until maturity). We consider 
individuals of two species, a and b. Individuals of species a are modelled with an 
higher growth rate with respect to individuals of species b (r a > rb). Carrying 
capacity for species a is, instead, lower than the carrying capacity for species 
b ( K a < Kb). The following CWC rules describe the r/K selection model for 
r a = 5 ,r b = 0.00125, K„ = 100 and K b = 1000: 


D 5 j 0.00125 L 7 

£ : a i —> a a £ : b \— > b b 

n 0.1 /ill 0.0000025 L 

£ : a a i —> a £ : b b \— > b 


We might consider a disruptive event occurring on average every 4000 years with 
the rule: 

T : (kJ A) i —> (x ] a b) 

devastating the whole content of the compartment (modelled with the variable 
A') and just leaving one individual of each species. In Figure [2] we show a 10000 
years simulation for an initial system containing just one individual for each 
species. Notice how individuals of species b are disadvantaged with respect to 
individuals of species a who reach the carrying capacity very soon. A curve 
showing the growth of individuals of species & in a stable (non disruptive) en¬ 
vironment is also shown. The full CWC model describing this example can be 
found at: http://www.di.unito.it/-troina/cmcl3/rK.cwc 
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Fig. 2. r/K selection in a disruptive environment. 


3.2 Competition and Mutualism 

In ecology, competition is a contest for resources between organisms: animals, 
e.g., compete for water supplies, food, mates, and other biological resources. 
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In the long term period, competition among individuals of the same species 
(intraspecific competition) and among individuals of different species ( interspe¬ 
cific competition ) operates as a driving force of adaptation, and, eventually, by 
natural selection, of evolution. Competition, reducing the fitness of the individ¬ 
uals involved p 2 ] has a great potential in altering the structure of populations, 
communities and the evolution of interacting species. It results in the ultimate 
survival, and dominance, of the best suited variants of species: species less suited 
to compete for resources either adapt or die out. We already depicted a form of 
competition in the context of the logistic model, where individuals of the same 
species compete for vital space (limited by the carrying capacity K). 

Quite an apposite force is mutualism, contest in which organisms of differ¬ 
ent species biologically interact in a relationship where each of the individuals 
involved obtain a fitness benefit. Similar interactions between individuals of the 
same species are known as co-operation. Mutualism belongs to the category of 
symbiotic relationships, including also commensalism (in which one species ben¬ 
efits and the other is neutral, i.e. has no harm nor benefits) and parasitism (in 
which one species benefits at the expense of the other). 

The general model for competition and mutualism between two species a and 
b is defined by the following equations l66j : 

"Y = r -ir^ ■ (K a -N a + a ab ■ N b ) 

^ • (K b - N b + a ba ■ N a ) 

where the r and I\ factors model the growth rates and the carrying capacities 
for the two species, and the a coefficients describe the nature of the relationship 
between the two species: if a,; ? - is negative, species Nj has negative effects on 
species N (i.e., by competing or preying it), if a t j is positive, species Nj has 
positive effects on species N (i.e., through some kind of mutualistic interaction). 

The logistic model, already discussed, is included in the differential equations 
above. Here we abstract away from it and just focus on the components which 
describe the effects of competition and mutualism we are now interested in. 

CWC Modelling 4 (Competition and Mutualism) For a compartment of 
type i, we can encode within CWC the model about competition and mutualism 
for individuals of two species a and b using the following stochastic rewrite rules: 

„ , fa-\aab\ f a a b if a ab > 0 „ , tvlaU f a b b if a ba > 0 

[b if a ab < 0 (a if a ba < 0 

where fi = is obtained from the usual growth rate and carrying capacity. The 
a coefficients are put in absolute value to compute the rate of the rule, their signs 
affect the right hand part of the rewrite rule. 

12 By fitness it is intended the ability of surviving and reproducing. A reduction in 
the fitness of an individual implies a reduction in the reproductive output. On the 
opposite side, a fitness benefit implies an improvement in the reproductive output. 
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Example 3. Mutualism has driven the evolution of much of the biological di¬ 
versity we see today, such as flower forms (important to attract mutualistic 
pollinators) and co-evolution between groups of species [57]. We consider two 
different species of pollinators, a and b, and two different species of angiosperms 
(flowering plants), c and d. The two pollinators compete between each other, and 
so do the angiosperms. Both species of pollinators have a mutualistic relation 
with both angiosperms, even if a slightly prefers c and b slightly prefers d. For 
each of the species involved we consider the rules for the logistic model and for 
each pair of species we consider the rules for competition and mutualism. The 
parameters used for this model are in Table [l] So, for example, the mutualistic 
relations between a and c are expressed by the following CWC rules 

‘&ac k c '&ca 

T : a c £—> a a c T : a c e—► a c c 

Figure [3] shows a simulation obtained starting from a system with 100 individ¬ 
uals of species a and b and 20 individuals of species c and d. Note the initially 
balanced competition between pollinators a and b. This random fluctuations are 
resolved by the “long run” competition between the angiosperms c and d: when d 
predominates over c it starts favouring the pollinator b that now can win its own 
competition with pollinator a. The model is completely symmetrical: in other 
runs, a faster casual predominance of a pollinator may lead the evolution of its 
preferred angiosperm. The CWC model describing this example can be found 
at: http://www.di.unito.it/~troina/cmcl3/compmutu.cwc 


Species (i) 

Ti 

Ki 

Ot-ai 

Oibi 

Ot-ci 

OLdi 

a 

0.2 

1000 

• 

-1 

+0.03 

+0.01 

b 

0.2 

1000 

-1 

• 

+0.01 

+0.03 

c 

0.0002 

200 

+0.25 

+0.1 

• 

-6 

d 

0.0002 

200 

+0.1 

+0.25 

-6 

• 


Table 1. Parameters for the model of competition and mutualism. 


3.3 Trophic Networks 

A food web is a network mapping different species according to their alimentary 
habits. The edges of the network, called trophic links , depict the feeding pathways 
(“who eats who”) in an ecological community [33]. At the base of the food 
web there are autotroph specie^] also called basal species. A food chain is a 
linear feeding pathway that links monophagous consumers (with only one exiting 
trophic link) from a top consumer, usually a larger predator, to a basal species. 

13 Self-feeding: able to produce complex organic compounds from simple inorganic 
molecules and light (by photosynthesis) or inorganic chemical reactions ( cheniosyn- 
thesis). 
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Fig. 3. Competition and Mutualism. 


The length of a chain is given by the number of links between the top consumer 
and the base of the web. The influence that the elements of a food web have 
on each other determine important features of an ecosystem like the presence 
of strong inter actors (or keystone species ), the total number of species, and the 
structure, functionality and stability of the ecological community. 

To model quantitatively a trophic link between species a and b (i.e., a par¬ 
ticular kind of competition) we might use Lotka-Volterra equations [71] : 

^=N b -(r b ~a- N a ) 

^=N a -(P-N b -d) 

where N a and N b are the numbers of predators and preys, respectively, r b is the 
rate for prey growth, a is the prey mortality rate for per-capita predation, jd 
models the efficiency of conversion from prey to predator and d is the mortality 
rate for predators. 


CWC Modelling 5 (Trophic Links) Within a compartment of type £, given 
a predation mortality a and conversion from prey to predator ft, we can encode 
in CWC a trophic link between individuals of species a (predator) and b (prey) 
by the following rules: 


£ : a 
£ : a 



b 


Here we omitted the rules for the prey exponential growth (absent predators) 
and predators exponential death (absent preys). These factors are present in the 
Lotka-Volterra model between two species, but could be substituted by the effects 
of other trophic links within the food web. In a more general scenario, a trophic 
link between species a and b could be expressed condensing the two rules within 
the single rule: 

7 


£ : a b 


a a 
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with a rate 7 modelling both the prey mortality rate and the predator conversion 
factor. 


Example f. Trophic cascades occur when predators in a food web suppress the 
abundance of their prey, thus limiting the predation of the next lower trophic 
level. For example, an herbivore species could be considered in an intermediate 
trophic level between a basal species and an higher predator. Trophic cascades 
are important for understanding the effects of removing top predators from food 
webs, as humans have done in many ecosystems through hunting or fishing ac¬ 
tivities. We consider a three-level food chain between species a, b and c. The 
basal species a reproduces with the logistic model, the intermediate species b 
feeds on a, species c predates species b: 


£ : a 1 


„ 0.0002 

£ : a a 1 —> a 


„ , 0.0004 

£ : a b 1 —!> 0 b 


„ , 0.0008 
£ : b c 1 —> c c 


Individuals of species c die naturally, until an hunting species enters the ecosys¬ 
tem. At a rate lower than predation, b may also die naturally (absent predator). 
An atom h. may enter the ecosystem and start hunting individuals of species c: 

0 °- 5 ? 0 u 0 0 3 ~r , , 1 0.003 , . .t 07 , 

t : c 1 — > • £ : b 1 —>■ • T : h [x J X) 1 —»■ [x J A h) £ : h c 1 — > h 

Figure |4] shows a simulation for the initial term h (*J 1000 * a 100 * b 10 * c) e . 
When the hunting activity starts, by removing the top predator, a top-down 
cascade destroys the whole community. The CWC model describing this example 
can be found at: http://www.di.unito.it/-troina/cmc13/trophic.cwc. 



Fig. 4. A Throphic Cascade. 


4 An application: Croton wagneri and Climate Change 

The knowledge of the relationships between the various attributes that define 
natural ecosystems (composition, structure, biotic and abiotic interactions, etc.) 
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is crucial to improve our understanding of their functioning and dynamics. Much 
in the same way, this knowledge could allow us to evaluate and predict the 
ecological impacts of global change and to establish the appropriate measures 
to preserve the intrinsic characteristics of the ecological systems minimising the 
environmental impact. 

Dry ecosystems are characterised by the presence of discontinuous vegeta¬ 
tion that may reflect less than 60% of the available landscape. The main pattern 
in arid ecosystems is a vegetation mosaic composed of patches and clear sites. 
Dry ecosystems constitute an ideal model to analyse the relationships between 
the attributes of the ecosystem and its functioning at different levels of organ¬ 
isation. This becomes particularly important for the crucial need, for this kind 
of ecosystem, to deal with environmental problems like the loss of biodiversity 
and climate change m- Morphological and functional changes shown by each 
species of these ecosystems are generally very different. As a consequence, the 
continuously adapting structure of the community will depend on the plastic¬ 
ity of each species in response to the major environmental changes. In the last 
decades, changes in the morphology of a species have been investigated from a 
functional perspective. This kind of analysis becomes particularly relevant be¬ 
cause, for example, the existence of plant species in environments with extreme 
climatic conditions may depend not only on the ability of the plant to struc¬ 
turally specialise but also on its capacity to adapt metabolically. Another key 
factor of dry ecosystems is the possibility to study special patterns of the vegeta¬ 
tion distribution to understand whether the plant communities show competitive 
or mutualistic relations. In particular, in dry (infertile) ecosystems, when com¬ 
petitive relations predominate, the spatial distribution of the vegetation tends to 
be uniform (with a quite regular distance between the plants); when mutualistic 
relations predominate, the vegetation tends to form clusters of plants; if none 
of these relations is observed, the spatial distribution of the vegetation is more 
random m- 

The study site is located in a dry scrub in the south of Ecuador (03°58'29" 
S, 01 o 25'22" W) near the Catamayo Valley, with altitude ranging from 1400m 
to 1900m over the sea level. Floristically, in this site we can find typical species 
of xerophytic areas (about 107 different species and 41 botanical families). The 
seasonality of the area directly affects the species richness: about the 50% of the 
species reported in the study site emerge only in the rainy season. Most species 
are shrubs although there are at least 12 species of trees with widely scattered 
individuals, at least 50% of the species are herbs. The average temperature is 
20° C with an annual rainfall around 600 mm, the most of the precipitation 
occurs between December and March, the remaining months present a water 
deficit with values between 200 and 600 mm of rain annually. Generally, this 
area is composed by clay, rocky and sandy soils |1]. 

In [ST] about 1300 different species belonging to the dry ecosystems in North¬ 
west South America have been identified. For this study we will focus on data 
available for the species Croton wagneri Mull. Arg., belonging to the Euphor- 
biaceae family. This species, particularly widespread in tropical regions, can be 
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identified by the combination of latex, alternate simple leaves, a pair of glands at 
the apex of the petiole, and the presence of stipules. C. wagneri is the dominant 
endemic shrub in the dry scrub of Ecuador and has been listed as Near Threat¬ 
ened (NT) in the Red Book of Endemic Plants of Ecuador |BH]. This kind of 
shrub could be considered as a nurse specie^] and is particularly important for 
its ability to maintain the physical structure of the landscape and for its contri¬ 
bution to the functioning of the ecosystem (observing a marked mosaic pattern 
of patches having a relatively high biomass dispersed in a matrix of poor soil 
vegetation) [4I] . 

In the study area, 16 plots have been installed along four levels of altitude 
gradients (1400m, 1550m, 1700m and 1900m): two 30mx30m plots per gradient 
in plane terrain and two 30mx30m plots per gradient in a slope surface (with 
slope greater than 10°). The data collection survey consisted in enumerating all 
of the C. wagneri shrubs in the 16 plots: the spatial location of each individual 
was registered using a digital laser hypsometer. Additionally, plant heights were 
measured directly for each individual and the crown areas were calculated accord¬ 
ing to the method in [63]. Weather stations collect data about temperatures and 
rainfall for each altitude gradient. An extract of data collected from the field can 
be found at: http://www.di.unito.it/~troina/croton_data_extract.xlsx. 
This data show a morphological response of the shrub to two factors: temper¬ 
ature and terrain slope. A decrease of the plant height is observed at lower 
temperatures (corresponding to higher altitude gradients), or at higher slopes. 


4.1 The CWC model 

A simulation plot is modelled by a compartment with label P. Atoms g , repre¬ 
senting the plot gradient (one g for each metre of altitude over the level of the 
sea), describe an abiotic factor put in the compartment wrap. 

According to the temperature data collected by the weather stations we corre¬ 
late the mean temperatures in the different plots with their respective gradients. 
In the content of a simulation plot, atoms t , representing 1°C each, model its 
temperature. Remember that, in this case, the higher the gradient, the lower 
the temperature. Thus, we model a constant increase of temperature within the 
simulation plot compartment, controlled by the gradient elements g on its wrap: 

T : ( x \ X) p ^{x\t X) p T : (g x\ t X) p "‘W 4 (g x\X) p 

Atoms i are also contained within compartments of type P, representing the 
complementary angle of the plot’s slope (e.g., 90 * i for a plane plot or 66 *i for 
a 24° slope). 

We model C. wagneri as a CWC compartment with label c. Its observed 
trait, namely the plant height, is specified by atomic elements h (representing 
one mm each) on the compartment wrap. 
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A nurse plant is one with an established canopy, beneath which germination and 
survival are more likely due to increased shade, soil moisture, and nutrients. 
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To model the shrub heights distribution within a parcel, we consider the 
plant in two different states: a “young” and an “adult” state. Atomic elements 
y and a are exclusively, and uniquely, present within the plant compartment in 
such a way that the shrub height increases only when the shrub is in the young 
state (y in its content). The following rules describe (i) the passage of the plant 
from y to a state with a rate corresponding to a 1 year average value, and (ii) the 
growth of the plant, affected by temperature and slope, with a rate estimated 
to fit the field collected data: 

c : y 1 — > a P :t i (x\y X) c \— > t i (x h\ y X) c 

4.2 Simulation results 

Now we have a model to describe the distribution of C. wcigneri height using 
as parameters the plot’s gradient (n * g) and slope (to * i). Since we do not 
model explicitly interactions that might occur between C. wagneri individuals, 
we consider plots containing a single shrub. Carrying on multiple simulations, 
through the two phase model of the plant growth, after 1500 time units (here 
represented as days), we get a snapshot of the distribution of the shrubs heights 
within a parcel. The CWC model describing this application can be found at: 
http://www.di.unito.it/~troina/cmcl3/croton.cwc. 

Each of the graphs in Figure [5] is obtained by plotting the height deviation 
of 100 simulations with initial term (n* g\m*i (*J y) c ) p . The simulations in 
Figures [5] (a) and (c) reflect the conditions of real plots and the results give a 
good approximation of the real distribution of plant heights. Figures [5] (b) and 
(d) are produced considering an higher slope than the ones on the real plots 
from were the data has been collected. These simulation results can be used for 
further validation of the model by collecting data on new plots corresponding to 
the parameters of the simulation. 

If we already trust the validity of our model, we can remove the correla¬ 
tion between the gradient and the temperature, and directly express the latter. 
Predictions can thus be made about the shrub height at different temperatures, 
and how it could adapt to global climate change. Figure [6] shows two possible 
distributions of the shrub height at lower temperatures (given it will actually 
survive these more extreme conditions and follow the same trend). 

5 Conclusions and Related Works 

The long-term goal of Computational Ecology is the development of methods to 
predict the response of ecosystems to changes in their physical, chemical and bio¬ 
logical components. Computational models, and their ability to understand and 
predict the biological world, could be used to express the mechanisms governing 
the structure and function of natural populations, communities, and ecosystems. 
Until recent times, there was insufficient computational power to run stochas¬ 
tic, individually-based, spatially explicit models. Today, however, some of these 
techniques could be investigated [53] . 


20 


P. Ramon and A. Troina 




(a) 1400 * g and 90 * i 


(b) 1550 * g and 60 * * 




(c) 1700 * g and 85 * i (d) 1900 * g and 75 * i 


Fig. 5. Deviation of the height of Croton wagneri for 100 simulations. 




(a) 12°C, plain terrain (b) 10°C, plain terrain 


Fig. 6. Deviation of the height of Croton wagneri for 100 simulations. 
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The main objective of Computational Ecology is twofold. First, we aim at 
discovering new models and theories within Computer Science inspired by the 
natural world, and at producing techniques and tools to deal with much more 
complex problems than those addressable with current technology. Second, we 
plan to provide ecologists with an environment for attacking problems at a sys¬ 
tem level, not addressable without using Information Technology. This envi¬ 
ronment will provide ecologists with modelling, analysis and simulation tools 
capable of handling complex behaviour and of representing emerging properties. 

Calculi developed to describe process interactions in a compartmentalised 
setting are well suited for the description and analysis of the evolution of eco¬ 
logical systems. The topology of the ecosystem can be directly encoded within 
the nested structure of the compartments. These calculi can be used to repre¬ 
sent structured natural processes in a greater detail, when compared to purely 
numerical analysis. As an example, food webs can give rise to combinatorial 
interactions resulting in the formation of complex systems with emergent prop¬ 
erties (as signalling pathways do in cellular biology), and, in some cases, giving 
rise to chaotic behaviour. 

Also, ecological niches describe how organisms or populations respond to the 
distribution of resources and competitors [451, and arise as the spatial sectors 
where organisms and populations tend to distribute, often forming geographical 
clusters of dominant species. Using the Spatial Calculus of Wrapped Compart¬ 
ments mm , we might be able to extend our analysis to the case of ecological 
niches and show how they could emerge from complex ecological interactions 
and then become a fundamental characteristic of an ecosystem. 

As a final remark about ecological modelling with a framework based on 
stochastic rewrite rules, we underline an important compositional feature. How 
can we test an hypothetical scenario where a grazing species is introduced in the 
model of our case study? A possibility could be to represent the grazing species 
with a new CWC atom (e.g. s) and then just add the new competitive rules to 
the previously validated model (e.g. the rule P : s (h x J X) c \—> s (a; J X) c ). 
Changing in the same sense a model based on ordinary differential equations 
would, instead, result in a complete new model were all previous equations should 
be rewritten. 

5.1 Related works 

Computer simulations play a rather unique role in Ecology, compared with other 
sciences, and there are several reasons for the extensive growth in the number 
of simulation -based studies in ecological mathematical modelling. In particular, 
the problem is that, although a field experiment is a common research approach 
in ecology, replicated experiments under controlled conditions (which is a corner¬ 
stone of all natural sciences providing information for theory development and 
validation) are rarely possible because of the transient nature of the environment: 
just consider, for instance, the impossibility to reproduce the same weather pat¬ 
tern for a repeated experiment in an ecological study. We also mention it here 
that large-scale ecological experiments are costly and, in the situation when 
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consequences are poorly understood, can have adverse effects on some species, 
put in danger the biodiversity of the ecosystem, and may even pose a threat to 
human well -being. Capturing the complexity of real systems through tractable 
experiments is therefore logistically not feasible. Mathematical modelling and 
computer simulations create a convenient “virtual environment” and hence can 
provide a valuable supplement, or sometimes even an alternative, to the field 
experiment. 

Ecological Modelling. It has long been recognised that numerical modelling 
and computer simulations can be used as a powerful research tool to understand, 
and sometimes to predict, the tendencies and peculiarities in the dynamics of 
populations and ecosystems. It has been, however, much less appreciated that 
the context of modelling and simulations in ecology is essentially different from 
those that normally exist in other natural sciences [ 551 - 

Ecology became a quantitative and theory-based science since the seminal 
studies by Lotka [33], Volterra m and Gause [3B], who were the first to use 
mathematical tools for ecological problems. General principles of ecosystem or¬ 
ganisation were later refined and systematised by Odum & Odum [49] . while 
the mathematical theory was further developed by Skcllam [64] and Turing [68], 
who emphasised the importance of the spatial aspect. The bright mathematical 
ideas of those seminal works sparked a huge fire. The last quarter of the twenti¬ 
eth century saw an outbreak of interest in mathematical ecology and ecological 
modelling [73151129146102100] . Especially over the last decade, more and more 
complicated models have been developed with a generic target to take into ac¬ 
count the ecological interactions in much detail and hence to provide an accurate 
description of ecosystems dynamics. Owing to their increased complexity, many 
of the models had to be solved numerically, a development that was inspired and 
made possible by the simultaneous advances in computer science and technology. 

The modelling approaches can be very different in terms of the mathematics 
used and depending on the goals of the study, and there are several ways to 
classify them. For instance, there is an apparent difference between statistical 
models [28] and “mechanistic” models [47150] . although simulation-based stud¬ 
ies may sometimes include both. Taken from another angle, two qualitatively 
different modelling streams are rule-based approaches (such as individual-based 
models and cellular automata) [33] and equation-based ones. 

Another way to sort out the models used in ecology is to consider the level 
of complexity involved. Depending on the purposes of the ecological study, there 
have been two different streams in model building [45j . In case the purpose is 
to predict the system’s state (with a certain reasonable accuracy), the model is 
expected to include as many details as possible. This approach is often called 
predictive modelling. The mathematical models arising in this way can be very 
complicated and analytically intractable in an exact way (in these cases the 
model can still be partially analysed via a limited number of runs of computer 
simulations) [53148154159] , Alternatively, the purpose of the study can be to 
understand the current features of the ecosystem, e.g. to identify the factors re- 
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sponsible for a population decline or a population outbreak, but not necessarily 
to predict their development quantitatively. We will call this approach a concep¬ 
tual modelling. In this case, the corresponding models can be pretty simple, even 
if their exact solutions are still not always possible; therefore, they often have 
to be solved by simulations as well I51I2QI62I . These two streams of theoretical- 
ecological research can be clearly seen in the literature, even though it is not 
always straightforward to distinguish between them as sometimes simple models 
may show a certain predictive power and, on the contrary, complicated ones are 
used for making a qualitative insight into some subtle issues. 

Formal Computational Frameworks. As P-Systems |58l59j and the Calculus 
of Looping Sequences (CLS, for short) [TT] . the Calculus of Wrapped Compart¬ 
ments is a framework modelling topological compartmentalisation inspired by 
biological membranes, and with a semantics given in terms of rewrite rules. 

CWC has been developed as a simplification of CLS, focusing on stochastic 
multiset rewriting. The main difference between CWC and CLS consists in the 
exclusion of the sequence operator, that constructs ordered strings out of the 
atomic elements of the calculus. While the two calculi keep the same expres¬ 
siveness, some differences arise on the way systems are described. On the one 
hand, the Calculus of Looping Sequences allows to define ordered sequences in 
a more succinct way (for examples when describing sequences of genes in DNA 
or sequences of amino acids in proteins) [^] On the other hand, CWC reflects 
in a more realistic way the fluid mosaic model of the lipid bilayer (for example 
in the case of cellular membrane description, where proteins are free to float), 
and, the addition of compartment labels allows to characterise the properties 
peculiar to given classes of compartments. Ultimately, focusing on multisets and 
avoiding to deal explicitly with ordered sequences (and, thus, variables for se¬ 
quences) strongly simplifies the pattern matching procedure in the development 
of a simulation tool. 

The Calculus of Looping Sequences has been extended with type systems 
in ltil.31132181161 . As an application to ecology, stochastic CLS (see H) is used 
in |12j to model population dynamics. 

P-Systems have been proposed as a computational model inspired by biolog¬ 
ical structures. They are defined as a nesting of membranes in which multisets of 
objects can react according to pre defined rewrite rules. Maximal-parallelism is 
the key feature of P-Systems: at each evolution step all rewrite rules, in all mem¬ 
branes, are applied as many times as possible. Such a feature makes P-Systems 
a very powerful computational model and a versatile instrument to evaluate 
expressiveness of languages. However, it is not practical to describe stochastic 
systems with a maximally-parallel evolution: exact stochastic simulations based 
on race conditions model systems evolutions as a sequence of successive steps, 
each of which with a particular duration modelled by an exponential probability 
distribution. 
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An ordered sequence can be expressed in CWC as a series of nested compartments, 
ordered from the outermost compartment to the innermost one. 
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There is a large body of literature about applications of P-Systems to ecolog¬ 
ical modelling. In [20121122] , P-Systems are enriched with a probabilistic seman¬ 
tics to model different ecological systems in the Catalan Pyrenees. Rules could 
still be applied in a parallel fashion since reduction durations are not explic¬ 
itly taken into account. In [13114115] . P-Systems are enriched with a stochastic 
semantics and used to model metapopulation dynamics. The addition of mute 
rules allows to keep a form of parallelism reducing the maximal consumption of 
objects. 

While all these calculi allow to manage systems topology through nesting 
and compartmentalisation, explicit spatial models are able to depict more precise 
localities and ecological niches , describing how organisms or populations respond 
to the distribution of resources and competitors [43]. The spatial extensions of 
CWC [T7j, CLS [5j and P-Systems [TO] could be used to express this kind of 
analysis allowing to deal with spatial coordinates. 
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